A Non-Demolition Single Spin Meter 
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We present the theory of a single spin meter consisting of a quantum dot in a magnetic field under microwave 
irradiation combined with a charge counter. We show that when a current is passed through the dot, a change 
in the average occupation number occurs if the microwaves are resonant with the on-dot Zeeman splitting. The 
width of the resonant change is given by the microwave induced Rabi frequency, making the quantum dot a 
sensitive probe of the local magnetic field and enabling the detection of the state of a nearby spin. If the dot-spin 
and the nearby spin have different g-factors a non-demolition readout of the spin state can be achieved. The 
conditions for a reliable spin readout are found. 

PACS numbers: 85.75.-d, 73.63.-b, 85.35.Gv, 75.50.Xx 



Single-shot readout of individual spins lies at the heart of 
many future technologies, from spintronics and quantum com- 
putation to single molecule spin resonance. There have been 
several recent proposals and demonstrations of single spin 
read out, all with restricted applicability, relying on the de- 
tails of the physical system in which the spin is housed fljjlor 
requiring certain special features, such as optical activity 12], 
nuclear spins H] or a large detector-system interaction Qj; 
others are too slow to achieve a measurement in a single at- 
tempt [5]]. Here we introduce a versatile spin meter that can 
achieve reliable determination of the spin state. Our electrical 
method is non-invasive, in contrast to standard techniques that 
destroy the spin in the course of the measurement [6]. It con- 
sists of a quantum dot under microwave illumination and can 
in principle be used to detect arbitrarily weakly coupled spins 
in one shot. 

Our scheme could be used to read out an arbitrary spin, and 
we will focus our discussion on a system that cannot be mea- 
sured by any means proposed so far: a molecular spin that is 
not optically active. Single molecule magnets exhibit a range 
of quantum phenomena 071] and represent a possible imple- 
mentation of spin qubits |8]. One of the leading candidates is 
a fullerene with an endohedral nitrogen (N@C6o), whose elec- 
tron spin has a remarkably long coherence time J9[]. However, 
the very property that makes the spin coherence time so long, 
i.e. weak interaction of the spin with the environment, also 
makes the readout difficult. One path to single spin readout 
is to place the molecule in a tiny gap between metal contacts; 
recent electrical current measurements through an endohedral 
fullerene show spin dependent features in the current voltage 
characteristic of the device fioll . The problem with this ap- 
proach is the unpredictable tunnel coupling of the molecule 
to the metal leads. In this Letter, we shall show that this dif- 
ficulty can be overcome by spatially separating the electrical 
readout from the molecular spin that is measured. 

The set-up for our proposed spin meter is shown schemati- 
cally in Fig. Q] A general obstacle in the sort of device we are 
proposing is the wide range of relevant energy scales, which 
at first seem to preclude sensitive spin measurement. An esti- 
mate of the dipole-dipole interaction energy between two elec- 
tron spins is J w 100 r/r o MHz, where ro = 1 nm. For a 
successful measurement, the current flow through the dot has 
to be sensitive to a corresponding energy change. The en- 
ergy resolution of the dot is usually governed by the width of 



Figure 1 : An implementation of the spin meter concept. The central 
part of the device is a quantum dot that is defined by electrical gates 
which pinch off a one-dimensional conductor, such as a nanowire 
or a carbon nanotube. The electrical transport through the quantum 
dot can be controlled with high precision by varying the source and 
drain potentials and the potentials of the extra finger gates placed at 
the junctions. The charging state of the dot can be sensed through 
an electrometer charge counter which is coupled to the quantum dot 
by a conducting bridge, and a back gate can control the number of 
electrons on the quantum dot. A single spin, e. g. in a fullerene, 
is placed close to the quantum dot. If there is no electron hopping 
between the dot and the fullerene, then the main coupling mecha- 
nism between the electron spin on the dot and the fullerene spin is 
the dipole-dipole interaction; there may also be an exchange interac- 
tion that could be tuned by inserting a conducting molecule between 
the fullerene and the nanotube. In such a setup no current will flow 
directly through the fullerene and thus the spin will not be disturbed. 



a typical electronic energy level, which depends on the tun- 
nelling rate and, more importantly, by the temperature depen- 
dent width of the Fermi function for the electrons in the leads. 
With attainable temperatures in standard dilution refrigerators 
(~ 10 mK) the energy resolution of the dot becomes ~ 1 
GHz, clearly not enough to resolve the dipole-dipole interac- 
tion. 

This obstacle can be overcome by considering spin reso- 
nance effects that are introduced by applying microwave ra- 
diation. Recently there has been an increasing interest in 
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methods for electrical detection of spin resonance ifTH [l2il . 
Mozyrksy and Martin showed that the current in a quantum 
channel next to an impurity spin depends on frequency of the 
incident microwave frequency [13]. We showed that a similar 
resonant change in the current takes place in transport through 
a quantum dot [14]. The line-width of the resonance intro- 
duces a new energy scale into the energy dependence of the 
dot occupations that can be much smaller than the temper- 
ature, subsequently enabling the detection of much smaller 
shifts in energy. Typical microwave intensities correspond to 
frequencies ~ 1 MHz but can be made smaller, provided that 
the lifetime of the electron on the dot is sufficient to allow for 
equalisation of the spin populations on resonance, i. e. the 
tunnelling rate should be smaller than the microwave induced 
Rabi frequency. Given a controllable tunnelling rate the sensi- 
tivity of the dot detector can be increased by lowering the mi- 
crowave intensity and also the tunnelling rate, with the lower 
limit set by the coherence time of the spin on the dot. 

The simplest model capturing all the relevant physics con- 
sists of two spin 1/2 particles, one situated on the quantum 
dot (the flying spin), the other nearby (the stationary spin). In 
a magnetic field the energy levels for spin up and spin down 
for each of the spins will be split by the Zeeman energy, but 
in general the splitting will be different for the flying and the 
stationary spin whose g-factors are typically unequal. Thus 
any incident narrow-band microwave radiation can only be 
resonant with one of the spins. We choose it to be resonant 
or nearly resonant with the flying spin. The dipole-dipole and 
exchange coupling between the two spins give rise to an effec- 
tive Ising interaction, since flip-flop processes are suppressed 
by the Zeeman terms. The resulting Hamiltonian in the rotat- 
ing frame after the rotating wave approximation can be written 
as 



H 



where S labels the stationary spin and F the flying spin, 
Aqs/f = Ss/fI^bBz — huio, with B z the magnetic field in 
the z-direction, hb the Bohr magneton, gs/F the g-factor for 
the stationary/flying spin respectively and ujq the microwave 
frequency. The Rabi frequency due to an oscillating magnetic 
field in the x-direction for the flying spin is SI; we can ne- 
glect the equivalent term for the off -resonant static spin. J is 
the Ising interaction strength. The Hamiltonian separates into 
two decoupled subspaces, associated with stationary spin up 
(|) and spin down (j) respectively, and is written as follows: 
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where we defined the spin subspace dependent detuning 
A j/j = Aqf T J- The upper/lower sign refers to the 
up/down subspace. 

The resulting energy level structure is shown in Fig.|2] De- 
creasing the magnetic field will shift the upper pair of lev- 
els towards resonance, with the Rabi frequency dominating 
their splitting. Meanwhile, the splitting of the lower pair in- 
creases. Close to resonance the eigenstates become products 
of the static spin state with linear superpositions of flying spin 
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Figure 2: The level structure of the two spin system in the rotating 
frame for Aof = 0. After decreasing the magnetic field the upper 
two levels (stationary spin down subspace) move into resonance. The 
difference in the Zeeman splittings is given by SA Z = Aos — Aof- 



up and flying spin down: 

l«>F = fl sT IT) ir + a A U)f 

\b) F = MT) F + Mi) F , 

s =|, I denotes the stationary spin subspace. The detuning 
dependent coefficients are given by 
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where fe s f = — a s j and b s i = a s -\. 

Having established a model we can now ask whether we 
retrieve information about the interacting two spin system by 
passing a current through the quantum dot while monitoring 
the average population on the dot. Considering only the sta- 
tionary spin up subspace (or similarly the spin down subspace) 
we can analyse the rates of tunnelling from the leads onto the 
dot and from the dot into the leads. A schematic depiction of 
the tunnelling rates is shown in Fig. [3] 

We can derive the rate equations for the stationary spin up 
subspace, as 



~ (l a l| 2 7| > + ktlSf) PTa, (1) 

where p^ a is the probability of finding the stationary spin in 
state | and the flying spin in state \a) f andpfo is the probabil- 
ity of finding the stationary spin in state j and no electron on 
the dot, jy^ are the tunnelling rates onto the dot from either 

spin up or spin down leads and 7^ are the tunnelling rates 
off the dot for different spin populations. They are given by 

7 TA = 7 ^T A + 7 «T A ' 7 fA = 7 ^T A + 7 «T A 

where the tunnelling rates from the left and right leads are 
given by 

it] = 7o/(>l - Thjo)) 7lt = 7o [1 - /Ol - fiwo)] 
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with similar equations for the right lead, /(e) = 1/(1 + 
exp(/3e)) being the Fermi function and the chemical potentials 
in the left and right lead are given by fj, L i R = fj, — V g ± V s d/2 
with /i being the chemical potential, V g the gate voltage and 
V s d the bias voltage. A similar rate equation holds for the 
population of other energy eigenstate 



(2) 



Additionally the individual probabilities have to add up to 
the probability to find the stationary spin in a spin up state 
P] = P]a + Ptb + P]o- The equivalent rate equations can be 
derived for the stationary spin down subspace and pf +p j = 1 . 
The stationary state can be obtained for each stationary spin 
subspace from pwb/o = and Pi a /b/o — and we obtain 
the average population on the quantum dot as a function of 
the applied magnetic field, 



l- PT L(A 0T )-p iJ L(A 
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The coefficients in the Lorentzian Eq. (0]l are given by 
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The behaviour of the average dot population, Eq. 0, as a 
function of the detuning is shown in Fig. |4] If we want to 
distinguish spin up from spin down we can monitor the popu- 
lation at resonance, say Aof = J and consider the change in 
population from a pure spin up state (p^ = 1) to a pure spin 
down state (pi = 1) 

An = L(0) - L(2J). 

The dependence of An on the spin-spin interaction strength 
J and the temperature is shown in Fig. [5] The change in occu- 
pation number decreases for increasing temperature and de- 
creasing spin-spin interaction strength. In order to detect a 
small change in the dot occupation number a certain num- 
ber of electron tunnelling events have to take place. To de- 
tect a change of An more than 1 / An 2 electron have to tun- 
nel, therefore the electron tunnelling rate needs to be much 
larger than the stationary spin relaxation time. In order to see 




Figure 3: A schematic diagram showing tunnelling off and onto the 
quantum dot. The states on the quantum dot are superpositions of 
spin up and spin down. The rates for tunnelling on and off the quan- 
tum dot depend on the spin, the bias voltage, temperature, Zeeman 
splitting and gate voltage. 




Figure 4: The average dot population (n) as a function of the detun- 
ing Aof for stationary spin up and spin down respectively. The two 
Lorentzian peaks are centred at ± J and have a width of a£l. Observ- 
ing the population at a detuning of + J the difference in population 
between the case where the stationary spin points up and where the 
stationary spin points down is given by An. 



a resonance of the spin on the quantum dot two things have 
to be fulfilled: The Rabi frequency has to be stronger than 
any spin relaxation, intrinsic or tunnelling induced, and the 
thermal spin polarisation has to be significantly different from 
zero. For a non-demolition measurement the difference in the 
spin Zeeman splittings between stationary and flying spin has 
to be much larger than the Rabi frequency and the spin-spin 
interaction. 

How does this translate to experimental parameters? With a 
spin coherence time of 1 /is on the carbon nanotube quantum 
dot we need a Rabi frequency of at least 17 = 10 MHz to drive 
the spin into saturation, and thus effect a change of the on 
dot spin occupation. Given the standard ESR microwave fre- 
quency of u>o — 10 GHz and the corresponding magnetic field 
of 350 mT, at dilution refrigerator temperature (T = 50 mK) 
we can read out a spin with a relaxation time of 0.1 ms and a 
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coupling of J = 3 MHz. To achieve such a coupling strength 
purely by dipole-dipole interaction the spin must not be fur- 
ther than 1 nm from the nanotube and the length of the dot 
must be « 100 nm. Such a setup would for example allow the 
readout of a molecular magnet spin [8J]. Using a conductive 
molecule to attach the spin to the nanotube could increase the 
coupling strength or allow larger distances between nanotube 
and spin. 

Quantum dots in nanotubes have been fabricated (see e.g. 
lfl5[ [IH]), also with integrated charge counters lfl7l [Hil . and 
on-chip microwave resonators reaching Rabi frequencies of 
~ 10 MHz have been demonstrated (e.g. EHIM]). We there- 
fore conclude that our scheme is realisable with current tech- 
nology. 



Figure 5: Change of the charge on the quantum dot between spin up 
and spin down on the stationary spin as a function of the exchange 
coupling between the stationary and the flying spin for small bias 
voltages (Vsd = O.lfl, V g = 0). Different curves are for different 
temperatures. The blue shading in the background indicates the de- 
tectable change in occupation given a certain relaxation time of the 
stationary spin. The green dot indicates the example discussed in the 
text. 
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